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ABSTRACT 


In this work, a physics-based simulation of non-ionizing proton radiation damage effects at different 
energy levels on a GaN-on-silicon high electron mobility transistor was created. Based on physical 
results of 2.0-MeV protons irradiation to fluence levels of 6 x lO^'^ protons/cm^, the simulation was 
tuned to match electron mobility and then compared to threshold voltage Vth on state resistance 
R„n and transconductance gm- A Monte Carlo simulator was used to model two particle interactions 
utilizing the Kinchin and Pease model. The model was developed in Silvaco ATLAS, but the 
Athena and Victory Stress modules were also utilized. After comparison of changing characteristics 
between the model and the physical device at 2.0-MeV proton irradiation, predictions were made for 
5.0, 10.0, 20.0 and 40.0-MeV proton irradiation. The model generally overpredicted damage in the 
lattice when compared to the physical results seen in prior work. 
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CHAPTER 1: 
Introduction 


The family of Ill-Nitride (III-N) semiconductor materials is used extensively in electronic 
devices requiring high-power density, high-temperature and high-frequency operation. Due 
to the wide band-gap and high mobility of III-N materials, devices fabricated using III-N 
materials outperform devices fabricated from related III-V semiconductors in breakdown 
voltage, radio-frequency (RF) power gain and thermal performance. The ability to fab¬ 
ricate III-N heterojunctions which lead to the formation of a high mobility, high density 
two-dimensional electron gas (2DEG) allow for the fabrication of high power RF transistors 
for use in applications where vacuum tubes have traditionally been the only option. Ad¬ 
ditionally, due to the high atomic threshold displacement energy of III-Ns, these materials 
promise superior reliability in environments where radiation effects are a concern. 

The 2DEG formed in III-N heterostructures can be leveraged to create a three-terminal 
field-effect transistor, known as a high-electron mobility transistor (HEMT), which is opti¬ 
mally suited for RE applications. III-N HEMTs, particularly those fabricated using (0001) 
oriented Gallium Nitride (GaN)-Aluminum Gallium Nitride (AlGaN) heterostructures, are 
currently found in many of our next generation weapons and radar systems. High bandwidth 
shipboard communications, phased-array radar and electronic warfare modules are exam¬ 
ples of systems that are already taking advantage of the properties wide band-gap III-Ns 
provide [I]-[3]. Due to the increased power density seen in GaN HEMTs, they are effective 
replacements for klystrons, which take significantly more space in shipboard applications. 
Another application for this technology is satellites. Again, because of their decreased size, 
modules leveraging GaN HEMTs save payload space without sacrificing capability. This is 
crucial for payload launch cost in space systems. The effects of radiation upon electronics 
in the space environment has been well documented over the years and provides engineers 
with different challenges as compared with terrestrial environments. Knowing how radi¬ 
ation affects device characteristics and reliability is of prime importance to the design of 
devices that will be used in space-based applications. Extensive research upon the effects 
of radiation upon the III-N material family has been carried out, but the understanding of 
the effects in devices, especially in GaN HEMTs, is not complete. 
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1.1 Research Objective 

The effects of non-ionizing damage induced by high dose massive particle radiation on GaN- 
on-silicon (GaN-on-Si) ffEMTs are explored using a physics-based technology computer- 
aided design (TCAD) model. The model is based upon GaN-on-Si HEMTs fabricated at the 
U.S. Naval Research Eaboratory (NRE) for radiation experimentation. The Silvaco ATEAS 
TCAD package is used to first simulate the transport of charge carriers through the initial 
GaN-on-Si HEMT structure based upon measurements taken at the Naval Postgraduate 
School (NPS). Special attention is paid to accurate modeling of the mobility of the 2DEG 
channel, which is degraded by the presence of radiation damage. 

Utilizing the stopping and range of ions in matter (SRIM) software package, we used a 
Monte Carlo simulator of heavy ion radiation damage within materials to generate damage 
profiles for interstitials and vacancies induced by proton irradiation of varying energies for 
the HEMT structure. The SRIM program utilizes the Kinchin and Pease (K-P) model of 
Erenkel pair generation to simulate damage events under heavy ion irradiation. Profiles 
were created for proton energies of 2.0, 5.0, 10.0, 20.0 and 40.0-MeV protons. These 
profiles were then inserted into the ATEAS model developed for the initial device. A 
literature review was conducted to identify and classify the energy levels of the related 
traps associated with vacancies and interstitials created by the radiation damage predicted 
by SRIM. Again, care was taken to select a physically realistic model of the mobility term 
which is affected by the presence of radiation induced defects. 

Results were compared with data taken from prior work at NPS and NRE involving 2.0- 
MeV proton irradiation of the NRE devices to the model outputs. After modeling the effects 
of 2.0-MeV proton irradiation, we expanded the model to predict damage characteristics at 
the additional proton energy levels, with emphasis on predicting the expected change from 
40.0-MeV proton irradiation for fiuence levels above 10^^ p’*' • cm“^. 


1.2 Related Work 

Studies conducted on proton radiation in AlGaN/GaN HEMTs often have a wide variety 
of design variation between the devices. The manufacturing process, materials used and 
dimensions of the device all play a large role in the characteristics of the device. While it is 
impossible to find independent research with the identical devices utilized, Hu et al. have 
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done similar studies on 1.8-MeV proton irradiation of AlGaN/GaN-on-sapphire devices that 
are applicable to the study of the effects of proton radiation damage at similar and increased 
energy levels [4]-[6]. Modeling of radiation effects has also been studied. Hui-Chia Yu 
wrote an extensive paper on concepts of modeling diffusion of diluted species to varying 
complexities [7], [8]. 

At NPS, research on the radiation hardness of AlGaN/GaN HEMTs has been ongoing 
since the early 2000s. Initial research was carried out on commercial off-the-shelf devices. 
Thesis work at NPS started with properties research and failure analysis before moving 
to irradiation effects around six years ago [9]. More recently, research has been focused 
on the physical inspection of radiation damage in NRL fabricated GaN-on-Si HEMTs 
and the irradiation of these devices. In conjunction with NRE, work was carried out to 
irradiate and characterize GaN-on-Si devices using 2.0-MeV protons. Wade examined 
the micro-structure of damage created due to the irradiation by 2.0-MeV protons at a 
fluence of 6.0 x 10^^ protons/cm^ [10]. In that study, upon physical inspection of imagery 
provided by scanning transmission electron microscopy, voiding was noticed under the gate 
fingers of the device. This voiding was not originally present in the device. Along with 
this observation, energy dispersive spectroscopy analysis preformed by Wessel revealed 
diffusion of Nickel into the Gold layer [11]. This diffusion at layer boundaries was also 
evident at the AlGaN/GaN interface. East neutron irradiation experimentation was also 
carried out. In an experiment by lobst, NRE GaN-on-Si HEMTs were irradiated with 
gamma radiation to document the effects, followed by irradiation with fast neutrons (between 
1.0 and 20.0 MeV) at the McClellan Nuclear Research Center. The effects of the gamma 
radiation were filtered from the effects of the fast neutron irradiation for analysis [12]. 

1.3 Thesis Organization 

This thesis is laid out as follows. The device characterization, basic physics of proton 
interactions in a crystalline structure and principles of vacancy and interstitial creation and 
diffusion are reviewed in Chapter 2. The integration of physical properties into the Monte 
Carlo simulation software SRIM along with development of functions of damage related to 
depth based on simulation results are discussed in Chapter 3. The process for implementing 
these functions in the TCAD software Silvaco ATEAS along with the changes made to 
mobility modeling in order to more accurately encapsulate all of the effects contributing 
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to the changes in mobility exhibited in the device are also developed. In Chapter 4, the 
implementation of the radiation effects into the TCAD model is presented. The results of 
the model in comparison to physical results seen from 2.0-MeV proton radiation is analyzed, 
and predictions are made upon the damage induced by higher energy, higher fluence proton 
irradiation. Conclusions and opportunities for further work can be found in Chapter 5. 
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CHAPTER 2: 
Background 


2.1 Properties of Gallium Nitride 

Device quality III-N material has been available for device fabrication since the 1990s. For 
power and RF applications, the most important material in the III-N family is GaN. GaN is a 
wide band-gap semiconductor with a band-gap of approximately 3.4 eV at 300 K, dependent 
upon its crystal structure. The crystal structure of GaN can be either zincblende or wurtzite. 
The atomic structure and band diagram for wurtzite GaN are shown in Figures 2.1 and 2.2, 
respectfully. Wurtzite’s crystal structure is hexagonal and is defined by two different lattice 
constants a and c. The wurtzite crystal structure is different from other crystal structures 
in that it lacks inversion symmetry. Due to this lack of inversion symmetry, GaN (and 
other III-Ns) are pyroelectric: they possess a net spontaneous polarization Psp aligned 
antiparallel along the hexagonal c-axis of the crystal. Additionally, III-Ns are piezoelectric 
and modulate the polarization of the material dependent upon the strain in the lattice. 



f Energy 



Figure 2.1. Crystal Structure. Source: [13]. 


Figure 2.2. Band Diagram. Source: [14]. 
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GaN is not an easy material to grow in bulk at reasonable cost, so it is often grown on a 
substrate. The use of a specific substrate can alter the material properties of the GaN. When 
creating devices, GaN is not grown directly on the substrates because differences in lattice 
constants result in strains that cause large defects in the GaN. This mismatch increases the 
dislocation density within the active region of the device, resulting in diminished device 
performance. For that reason, buffer layers are grown on the substrate before growing GaN. 
Each buffer layer absorbs a portion of the stress so that when GaN is deposited, the stress 
is at an optimal level for crystalline growth. Traditional substrates for III-N growth, such as 
sapphire and silicon carbide, can be expensive materials to manufacture and do not have the 
technological maturity of older materials such as Si. The desirability of Si as a substrate 
for growth stems from the ability to grow cheap, large, high quality substrates; however, 
the lattice mismatch between Si and III-Ns is severe. This necessitates the use of several 
additional buffer regions for successful III-N growth. In GaN-on-Si devices studied in prior 
research at NFS, aluminum nitride (AIN) in combination with varying mole fractions of 
AlGaN were used as transitional layers before the GaN layer. 

The intrinsic strain developed in a crystal lattice is the result of an offset in the lattice 
constant between the growth layer and the substrate. This is referred to as lattice mismatch 
strain 

e=-, (2.1) 

at 

defined by the lattice constants of the substrate as and the lattice constant of the layer making 
contact with it ai- While too much strain can be a major issue during crystal growth, strain 
can be engineered to be useful in III-N structures. Once a high-quality III-N layer is grown, 
additional layers of III-Ns can be grown pseudomorphically, in which the entire epitaxial 
layer is under strain imposed via growth atop the substrate III-N layer. These additional 
layers create high quality heterojunctions which can be used for device applications. 

2.2 Spontaneous and Piezoelectric Polarization 

Polarization in III-Ns manifests as both spontaneous Psp and piezoelectric Pp^ polarization. 
The sum of these is the total polarization 


p - p + p 
^ sp ^ ^pz* 


( 2 . 2 ) 
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The units of polarization are C/cm^. The spontaneous polarization can be calculated 
for III-N alloys by linearly interpolating between the known values of the two composite 
materials of the alloy. For example, for AlGaN of mole fraction jc, 


Psp,AlGaN{^') — Psp,AlN^ "I” Psp,GaN{^ 


(2.3) 


Piezoelectric polarization is defined as 

Pp, = 2 |ei3 - e, (2.4) 

where the value eij comes from the piezoelectric coupling matrix and has units of C/cm^. 
The Cij terms come from the stiffness matrix of the strained material and have units of 
GPa[15]. 

In real heterostructures, dislocations due to imperfect growth of the III-N substrate layer 
relax the strain in subsequent layers. This leads to a reduction in piezoelectric polarization. 
The degree of relaxation can be quantified by the unitless coefficient rj, which varies between 
0 and 1. The relaxation coefficient can be used to calculate the reduction in polarization in 
the heterostructure due to strain relaxation [15] 


PT = Psp + {l-ri)Pp,. (2.5) 

The net effect of polarization is equivalent to the placement of a fixed sheet charge at the 
endpoints of the material layer. In a heterojunction, the net fixed charge is determined by the 
polarization difference between the two materials. This fixed charge is depicted in Figure 
2.3(b). The sign of the charge at the interface depends upon the growth direction of the 
heterostructure. In an AlGaN/GaN heterostructure grown in the (0001) direction, this sheet 
charge is positive. In order to maintain net charge neutrality, electrons accumulate to screen 
out this sheet charge at the interface. This is called a 2DEG due to the fact that the large 
electric field in the GaN results in the electrons being quantum mechanically confined in a 
triangular quantum well. The energy band diagram of such a structure is shown in Figure 
2.3, where the presence of the 2DEG can be clearly seen where the Eermi level goes above 
the conduction band in the GaN layer. 
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Conduction Band 
Fermi Level 

Valence Band 


(b) 


Figure 2.3. Electronic Structure of an AlGaN-GaN Heterostructure. Adapted 
from [15]. 


The band diagram and the sheet charge distribution of a representative AlGaN-GaN 
heterostructure are shown in Figures 2.3(a) and 2.3(b), respectively, where eN’^pp is the 
charge from surface defects at the top of the AlGaN layer, -Q^(AlGaN) is the polarization 
of the AlGaN at the top of the layer, -l-Q^(AlGaN) and -Q^(GaN) are the polarizations at 
the AlGaN/GaN interface, -l-Q^(GaN) is the polarization at the bottom of the GaN region, 
and -i-Qscr is the charge from defects at the GaN-substrate interface. 
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2.2.1 DC Characteristics of GaN HEMTs 

A HEMT is a three-terminal deviee based upon the modulation of the eonduetivity of a 
2DEG ehannel in a heterostrueture. Sueh a deviee is easily fabrieated through the deposition 
of a Sehottky barrier upon the heterostrueture and two ohmie eontaets to eontaet the 2DEG. 
In an ungated heterostrueture the 2DEG eannot be modulated direetly. The Sehottky gate 
eontaet is used to eontrol the eoneentration of the 2DEG and ean deplete the 2DEG given 
enough applied gate voltage. 


The DC eharaeteristies of a HEMT gatestaek ean be analyzed through energy and 
eharge balanee eonsiderations. Erom sueh an analysis, the values for the gate threshold 
voltage [16] 

-^Ec AP 

Vth = <l>b -^ ^ (2.6) 

<? '-eff 

neeessary to deplete the 2DEG ehannel ean be ealeulated, where (pb is the barrier height 
of the Sehottky eontaet, AEc is the eonduetion band energy level differenee between the 
AlGaN and GaN layers, AP is the differenee in polarization between the AlGaN and GaN 
layers, C'^ff is the eapaeitanee due to the thiekness of the AlGaN layer and the 2DEG offset 
due to quantum meehanieal eonfinement of the 2DEG, and q is the eleetron eharge. The 
equilibrium 2DEG eoneentration ean be ealeulated from [16] 


^2deg ~ 



q 


(2.7) 


The DC drain eurrent for GaN HEMTs 

/ y/ 2\ 

iv;,, - VnK, --f I ( 2 . 8 ) 

where W is the width of the deviee, L is the gate length, // is the mobility, Cgff is the 
eapaeitanee due to the thiekness of the AlGaN layer and the 2DEG offset due to quantum 
meehanieal eonfinement of the 2DEG, Vg^, and are the voltages at the points shown 
in Eigure 2.4, and Vth is the threshold voltage, whieh deseribes the drain eurrent dependenee 
upon the loeal voltages beneath the gate eontaet and does not take into aeeount resistanees 
that are present in the souree and drain aeeess region of the deviee. 
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Based on the diagram of the deviee shown in Figure 2.4, Equation 2.8 ean be modified 
based on voltage loop calculations to include the effects of the access resistances. 



Figure 2.4. Source and Drain Resistance Diagram 


The relationship between the applied drain and gate voltages V^s and Vgs and the 
gate-edge referenced drain and gate voltages and Vg^ are given by 


Vds = IRs + IRd + Vi, 


(2.9) 


and 

Vgs=IRs+V'g,. (2.10) 

From considerations of the 2DEG equilibrium concentration and electron mobility, the 
value of the resistances in the two access regions around the gate can be derived as [16] 


Rd 


^gd 


( 2 . 11 ) 


where Lgd is the distance between the gate edge and the drain edge. These parameters are 
important to the device and affect the characteristic curves of the device. The resistance 
values in these regions have a large affect on the device at higher drain currents. 


2.3 Non-Ionizing Radiation Damage Effects 

When an ionized particle moves through a crystal lattice, it interacts with the atoms in the 
lattice. These particles may directly collide with atomic nuclei physically or may "collide" 
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with atoms through charge interaction. The probability of a high energy ion colliding with 
atoms in the lattice is inversely proportional to the energy level of the ion [17]. When a 
collision does occur, conservation of energy and momentum is maintained and there can be 
a number of outcomes. If the particle has enough energy to overcome the energy required 
to displace the atom, known as the threshold displacement energy Eth, it knocks the lattice 
atom out of place. This creates an ionized atom that moves through the lattice with an 
energy dependent upon the nature of the collision, and a vacancy is left where it was before 
the collision. This atom is known as a primary knock-on atom (PKA). Depending on the 
amount of energy transferred to the PKA, we know it either comes to rest or continues to 
move. If it comes to rest, it either becomes an interstitial atom or replaces a vacancy in 
what is known as a replacement collision. If the atom continues to move, it may go on 
to collide with additional atoms. The atom continues to move and collide until its energy 
decreases below Eth, creating additional damage. Interstitial and vacancies remaining after 
the damage event are known as Frenkel pairs. A simple depiction of this process is seen in 
Figure 2.5. 
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Figure 2.5. Depiction of Particle Collision in a Crystal Lattice 
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The radiation damage events are classified as two-particle scattering events, involving 
only the incident ion and the impacted lattice atom. In lower energy collisions, there may 
be more than two particles involved in a collision event. As the majority of non-ionizing 
radiation damage of concern in semiconductors is induced by high energy particles, this 
work only considers the effects of two-particle interactions. 

2.4 Non-Ionizing Radiation Damage in GaN 

Exposing the semiconductor materials to non-ionizing radiation damage causes several 
different effects. The first is damage to the lattice. When atoms are displaced from their 
ordered lattice locations to a more disordered location, electron trap states are created. 
There are several types of trap states. They can be either acceptor or donor traps. Both 
types can potentially accept or donate more than one electron or hole. If a trap state is filled, 
it is said to be ionized. If enough trap states are ionized in a region, this results in charge 
build up, which can alter the field levels within a device and lead to shifts in the device 
characteristics. This can lead to undesired changes in circuit operation. 

In GaN layers, three trap states are known to be common after proton irradiation of 
fluence levels above 10^^ protons/cm^. These are referred to as ERl, ER2 and ER3. The 
"E" stands for electron traps, and the "R" indicates that it was developed from high-energy 
particle irradiation. Work by Auret in [18] has found that these three energy levels are found 
at 0.13,0.16 and 0.2 eV below the conduction band of GaN after 2.0-MeV proton irradiation. 
Auret does not specifically attribute these defects to specific defect types or atoms, and they 
are all shallow traps. In other work by von Bardeleben, a nitrogen split interstitial defect is 
noted in GaN after 12.0-MeV proton irradiation to fluence levels of 8 x 10^^ protons/cm^ 
which occurs around 1.0 eV below the conduction band by using electron paramagnetic 
resonance [19]. Most notably, a review of radiation effects by Pearton [20] examined and 
attempted to fingerprint individual radiation induced defects in the GaN lattice and the trap 
energy levels associated with them. In general, radiation induces nitrogen vacancies which 
act as shallow donors at 0.06 eV below the conduction band, nitrogen interstitials acting 
as deep acceptors around 1.0 eV below the conduction band, Ga vacancies forming deep 
acceptors at 1.0 eV above the valence band, and Ga interstitials forming deep donors at 0.8 
eV below the conduction band. 
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The crystal lattice is constantly vibrating and moving because of thermal energy 
maintained by the atoms and their bonds. This allows atoms in the lattice to change their 
position without any outside influence. If the temperature is not to the point of melting 
the material, the device can anneal out vacancies and interstitials. After irradiation events, 
devices try to return to their original form, and device characteristic have been seen to 
recover from these effects even at temperatures as low as 100 K [21]. In this work. He, Ar, 
O, Au, Te and Xe ions were implanted in GaN grown on sapphire at 15 K and Rutherford 
Backscattering analysis was performed to identify movement of ions as temperature was 
increased. 

Sometimes, when radiation displaces an atom, the atom is transferred enough energy to 
move it substantially from its vacancy. This makes it difficult for it to return and anneal out 
of existence or return to ordered form. When larger amounts of damage are created, piles 
of interstitials and vacancies start to form, and irreversible voiding occurs in the material. 
This has been seen in work at NFS in [10] and work by Limpijumnong in [22]. In [10], 
Wade visually inspected a HEMT that had undergone 2.0-MeV proton irradiation and found 
an accumulation of vacancies manifesting as a void under the Ni gate finger of the device. 
Limpijumnong’s work sought to classify and explain native defect diffusion in GaN crystal 
lattice through density functional theory. 
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CHAPTER 3: 
Model Development 


3.1 Device Composition 

The GaN-on-Si wafers for the HEMT used in this study were grown by Nitronex Inc. The 
substrate is 350-jum thick Si. A silicon nitride (SiNx) interface layer with a width of 2.0-nm 
and an AIN buffer layer of 300.0 nm were grown as an initial growth buffer layer, followed 
by two layers of AlGaN. The first is 400.0-nni thick with an A1 mole fraction of 0.7, 
while the second is 300.0 nm thick with an A1 mole fraction 0.4. A 0.8 /im high quality 
GaN layer was grown on top of the compliance layer. A 17.5-nm thick AlGaN layer with 
an A1 mole fraction of 0.27 was grown pseudomorphically as a barrier layer to form the 
2DEG channel. HEMT Device fabrication was carried out at NRE. The device gatestack 
consists of a 30.0-nm thick Ni gate, followed by a gold overlay of 300.0-nm. The device is 
passivated with a 100.0-nm thick SiNx layer. Ohmic contacts to the 2DEG were formed via 
the deposition of a Ti/Al/Ni/Au alloy followed by rapid thermal annealing. The diagram of 
the HEMT gate stack structure is shown in Eigure 3.1. 



Al. .Ga, ,N 


Si (111) Substrate 


AIq 27^30.73!' 


^13^4 


Eigure 3.1. Representation of the GaN High-Electron Mobility Transistor 
Gate Stack Structure 
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3.2 Initial Device Characterization 

The initial device electrical characterization was carried out using a probe station and a 
Keysight B1500A semiconductor parameter analyzer. DC IV measurements were per¬ 
formed to characterize the Ids versus Vgs transfer characteristics, Igs versus Vgs gate leakage, 
and Ids versus Vds characteristics of the device. Thirty-six devices on three identical reticles 
were tested to identify a device that could serve as a baseline for the development of the 
non-irradiated portion of the model. The section of the reticles considered had geometries 
identical to those used in 2.0-MeV irradiation experiments performed upon these devices 
in past work. The device dimension of the selected device are illustrated in Figure 3.2 in 
units of /rm. 



Figure 3.2. Scanning Transmission Electron Microscopy Photograph of a 
Single High-Electron Mobility Transistor from Reticle. Adapted from [11]. 


3.2.1 Initial Silvaco Model Development 

The initial model of the device was created in Silvaco ATEAS based on the physical 
dimensions of the actual device as seen in Eigure 3.3. 
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Figure 3.3. Device Created In ATLAS versus Scanning Transmission Electron 
Microscopy Imaging of Device. Adapted from [11]. 


The mesh was emphasized around the GaN and AIGaN interface and vertically under 
the edges of gate. Initially, a built-in polarization statement, polarization, was used 
for the AIGaN and GaN layers to simulate the spontaneous and piezoelectric polarization. 
Neither the AIN nor the AIGaN buffer layers received polarization statements because of 
the defective nature of these layers in the device. These layers are assumed not to provide 
any net polarization. Using these basic model parameters, we extracted the band diagram 
seen in Figure 3.4, showing the presence of a 2DEG at equilibrium. 
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AIGaN/GaN Band Diagram 



Figure 3.4. Silvaco Generated Band Diagram Beneath the Gate Staek 


The polarization statement was replaeed with a sheet eharge statement in the final 
version of the model. This provides the same effeet and does not alter the models eharae- 
teristies. The relaxation eoeffieient rj, defined in Equation 2.5, was adjusted to ealeulate the 
appropriate sheet eharge to fit 2DEG density measured via Hall measurements [23]. In or¬ 
der to model heating throughout the deviee, a temperature dependent thermal eonduetivity 
model was utilized, implemented via TCON. POWER from ATEAS’ GIGA self-heating mod¬ 
ule, for eaeh material in the deviee. The eharaeteristie equation for the thermal eonduetivity 
model is 


k{T) = 


a 

(7V3^‘ 


(3.1) 


The units for k are W/emK and the eoeffieients a and /3 for their respeetive layers are found 
in Table 3.1. 
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Table 3.1. Thermal Heating Coefficients 


Material 

a 

j8 


1.381 

0.317 

GaN 

0.65 

0.28 

Alo^GaojN 

0.688 

1.5535 

AIN 

1.425 

1.64 

Si 

1.48 

0 

Au 

2.0 

0 

Ni 

1.58 

0 


The source, drain, and gate were added to the simulation with electrode statements. 
The source and drain were modeled as ohmic contacts with work functions of 3.6 eV 
to match the affinity of the AlGaN and with an intrinsic 10.0-Q resistance. The 10.0-Q 
resistance was added to appropriately simulate the measured contact resistance of the ohmic 
contacts. The gate was modeled as a Schottky contact with a barrier height of 1.1 eV. The 
Shockley Reed Hall model was used to simulate carrier recombination in the model. 

3.2.2 Mobility Modeling 

Using appropriate mobility models is a key aspect in accurately modeling III-N devices. An 
initial, static value of 1260 cm^/(V • s) was adopted from Hall measurements taken on the 
bare reticle [23]. We quickly noted that, based on measured TV characteristics, this value 
resulted in a large overestimate of device currents. Comparison between measured device 
characteristics and the results of this initial mobility model are shown in Figures 3.5 and 
3.6. 
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Figure 3.5. Vg versus Id Vd = 100.0 mV 



Figure 3.6. Vg versus Id Vd = 10.0 V 
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Initial attempts at accurate mobility modeling were undertaken by implementing a 
custom mobility model within the GaN layer of the device. ATLAS allows users to create 
a custom mobility model which can be implemented via a c-function through the tofimun 
flag in the mobility statement. As a starting point, the well-known Albrecht model 
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(3.2) 


was utilized for GaN low-field mobility [24]. In the Albrecht model, ij.{N,Tl) is the 
mobility as a function of doping and lattice temperature, N is the total doping concentration 
in cm“^. No is the initial doping concentration in cm“^ and Tl is the lattice temperature in 
K. Parameters A, B and C are scaling factors to adjust the impact of the individual terms 
of Equation 3.2. It is important to note that Equation 3.2 was developed for bulk GaN. 
Because the Albrecht model was derived for bulk material, it is technically incorrect to use 
it to model electron mobility in the 2DEG. The first term of the model. 


A-N 





(3.3) 


accurately models impurity scattering in the bulk but not in the 2DEG. The second and third 
term, used to model acoustic and polar optical phonon scattering, can be shown to have the 
same temperature dependence as the corresponding mobility terms in the 2DEG-specific 
model [25]; therefore, those terms were incorporated as is and provided the foundation for 
the model. Charge dislocation and alloy scattering are additional factors that this model 
does not individually tailor terms for but can be shown as well to be constants independent 
of temperature and, thus, are incorporated in the tuning of the Albrecht model’s B and C 
terms. The values for the B and C coefficients were initially scaled to give an initial mobility 
matching that of the measurements in [26]. 


Eor high field mobility, the Barnes model was utilized [27]. The Barnes model is 


MnO+(Vsat/E)(E/Ecrity 
1 + (E/E,rity 


(3.4) 
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where yU„o is the low field mobility, Vsat is the saturation veloeity, E is the eleetrieal field, Ecrit 
is the critical electrical field value, and y is a unit less exponential fitting term. The Barnes 
model derives from measurements of the 2DEG electron velocity field scattering [27]. 
Values used for Vsat and Ecrit were 3.0 x 10^ and 20.0 x 10^, respectively. Based on recent 
work, other high field models showing velocity overshoot under conditions of high electric 
field may not be accurate [28]. In the work by Bajaj, high fields were utilized to try to reach 
velocity overshoot but found that electron velocity did not overshoot but merely saturated. 
More importantly, mobility overshoot was found to occur. This model also exhibits mobility 
overshoot, which fits the physical characteristics more accurately. 

Although modulation of the Albrecht scale factors brought the simulated output cur¬ 
rents into line with measured values, significant difficulty was found matching simulated 
to measured IV characteristics in regions of bias where the current was limited by access 
resistances as seen in Figures 3.7 and 3.8. To mediate this difficulty, measurements were 
taken to determine the access resistance for the gate and source regions of the device and 
can be seen in Figure 3.9. 



Figure 3.7. Vg versus Id Vd = 100.0 mV 


22 









Drain Current (A) 



Figure 3.8. Vg versus Id Vd = 10-0 V 



Gate Current (A) 

Figure 3.9. Source and Drain Resistance Measurements 
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From Equation 2.11, the width in the souree and drain regions along with the n 2 deg 
values should be the eonstant in a deviee. This leaves the mobility term and the length of 
the region that differs between the souree and drain. If the ratio of the resistanees is not 
equal to the ratio of the aeeess region lengths, then the mobility term must be the souree of 
variation. The ratio of resistanees in this deviee is 280/90, whieh is not equal to the ratio of 
the aeeess region lengths, 13/2. 

From this information, an additional faetor whieh eould explain the measured deviation 
of aeeess region mobility from expeeted values was sought. Reeent work by Lin to explain 
the diserepancy in [29], [30], proposed that the mobility in the aeeess regions is redueed due 
to an effeet known as polarized eoulomb field (PCF) seattering. PCF seattering is due to the 
faet that strain in the AlGaN-GaN interfaee region is not eonstant throughout the deviee. 
The strain is ehanged by the deposition of materials such as the source, gate and drain 
metal and the passivation layer. This alters the electron mobility as a function of distance 
between the source and drain of the device. As a result, the mobility model was expanded 
to include this effect. In order to properly account for PCF scattering potentially induced 
by deposition of additional layers above the AlGaN, the Silvaco model was shifted from 
ATLAS to Athena in order to utilize Silvaco’s Victory Stress functionality. After running 
the model in Victory Stress, non-uniform strain under the nitride layer was apparent as seen 
in Figure 3.10 which shows the XX component of strain in the lattice. The variation of 
the strain causes scattering and a decrease in electron mobility from the measured values 
predicted by Hall measurements on a bare heterojunction. 
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Figure 3.10. Athena Vietory Stress Simulation 


Although Vietory Stress eould not be used to direetly prediet the elFeets of inhomo¬ 
geneous strain in the AlGaN layer upon the mobility, the effeets of PCF seattering were 
implemented in another way. In order to take into aeeount the PCF seattering, the GaN 
region of the deviee was split into three seetions. This allowed us to eonstruet slightly 
different mobility terms for eaeh region based on the resistanees derived for eaeh region. 
The resulting eonstants used to fit the mobility are presented in Table 3.2 


Table 3.2. Constant Values for GaN Regions 



GaN Region 0 

GaN Region 1 

GaN Region 2 

A 

0.25 

0.25 

0.25 

B 

4.0 

1.35 

2.0 

C 

4.0 

1.35 

2.0 
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3.2.3 Final Model Evaluation 

Our initial comparison between the physieal deviee and the simulated deviee did not inelude 
any radiation effeets and was used as a measure for the aeeuraey of our starting model. After 
sealing the model, a strong eorrelation to the physieal deviee is shown in the eharaeteristie 
IV eurves illustrated in Figures 3.11, 3.12 and 3.13. 



Figure 3.11. Vg versus Id Comparison Between Model Output and Measured 
Results, Vd = 100.0 mV 
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Figure 3.12. Va versus Id Comparison Between Model Output and Measured 
Results, Vd = 10.0 V 



Figure 3.13. Id versus Vd Comparison Between Model Output and Measured 
Results, Vg = -3.0, -2.0, -1.0 and 0.0 V 
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3.3 SRIM Model 

In order to model the proper amount of radiation damage in the AlGaN/GaN-on-Si HEMT, 
the Monte Carlo ion radiation damage simulator SRIM was utilized . A layered device that 
has the equivalent dimensions of a two-dimensional slice of the HEMT underneath the gate 
was designed. These dimensions are found in Table 3.3. 

Table 3.3. Device Dimensions Eisted from Top to Bottom 


Material 

Thickness 

5 / 3 A 4 

100.0 nm 

Au Gate 

300.0 nm 

Ni Gate 

30.0 nm 

A/o.27Cao.73-^ 

17.5 nm 

GaN 

800.0 nm 

AIq/^Guq^N 

300.0 nm 

AIqjGqq^iN 

400.0 nm 

AIN 

300.0 nm 

SZ 3 A 4 

2.0 nm 

Sz’(lll) Substrate 

350.0 yum 


After the appropriate dimensions and material names were placed in SRIM, the correct 
densities were added. Eor layers that consisted of a monochromatic material, the default 
density was used. Eor layers that have specific molar compositions, the density was calcu¬ 
lated. The compound correlation value was left at one. The symbols, name, atomic number, 
weight, atom stoichiometry. Damage Eatt and Surf were all left at their default values. The 
displacement threshold energy values found by Pearton were used [31]. It is noted that while 
the displacement threshold energy levels for the same atom vary in different compounds, 
the assumption is made that in GaN, AlGaN, and AIN the thresholds were relatively similar 
and, therefore, were left as the same. The material properties that were changed can be 
found in Tables 3.4 and 3.5. The type of TRIM calculation used was a "Detailed Calculation 
with Eull Damage Cascades." The ions used to simulate p"*" irradiation were hydrogen atoms 
with the default mass and charge, zero angle of incidence, and 2.0-MeV energy level. 
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Table 3.4. TRIM Compound Setup Values 


Compound 

Density (g/cm^) 

5/3 V 4 

3.2 

Au 

19.311 

Ni 

8.8955 


3.0327 

GaN 

6.5 

AIq 4 GCIQ (jN 

2.8246 

AlojGao^^N 

2.3443 

AIN 

3.26 

5/3 V 4 

3.44 

Sz(lll) 

2.3212 


Densities of eomposite materials were linearly interpolated based on molar eomposi- 
tion. 


Table 3.5. TRIM Atomic Setup Values 


Atoms 

Displacement 

Threshold 

Energy (eV) 

Si 

15.0 

N 

28.0 

Au 

25.0 

Ni 

25.0 

Al 

72.0 

Ga 

72.0 

N 

32.4 

S/(lll) 

15.0 


Displacement threshold values from [31] 
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The SRIM program only collects data in 100 bins at a time regardless of the size of the 
device being modeled. The larger the structure, the wider the bin is required to be, reducing 
the resolution of the data collected. Due to the size variation of layers within the device, 
several simulations were run with windows of varying widths. Emphasis was placed on 
boundary regions and areas where current flows within the device. The following window 
sizes were utilized in our simulations: 0-3500, 3500-6500, 6500-22475 and 22475-3522500 
angstroms. The window size has no effect on the simulation calculations, it merely allows 
the bin size to be manipulated. 

Additional simulations were run at energy levels of 5.0 MeV, 10.0 MeV, 20.0 MeV and 
40.0 MeV. Due to the decreasing amount of particle interaction with the lattice as energy 
level increases, the number of ions simulated increased for increasing energy levels. In order 
to get accurate information about the interaction of 2.0-MeV proton radiation, 10 million 
particles were simulated. For 40.0-MeV protons, 50 million particles were simulated. 

The resultant damage through irradiation was estimated through measurements of 
interstitials and vacancies from the materials in each window. Energy levels less than 20.0 
MeV showed a linear trend in the interstitial and vacancy formation in both the AlGaN and 
GaN regions. The trend was directly correlated to the thickness of the layer. The linear fit 
for the AlGaN and GaN regions at 2.0-MeV protons are shown in Tables 3.6, 3.7 and the 
interstitial distribution in the AlGaN and GaN regions in Figures 3.14 and 3.15. The 20.0 
and 40.0-MeV energy levels showed a constant trend of damage in the lattice. The fitting 
information for the additional energy levels is found in Appendix, Section A. 1. 
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Interstitials per Ion per cm^ (Int/Ion/cm^) 



Figure 3.14. Linear Fitting for 2.0-MeV Protons in GaN. 


Table 3.6. Interstitial and Vacancy Distribution Fitting Functions in the GaN 
Region of the Device 


Lines for 2.0-MeV energy levels. 


Atom 

Slope 

F-intercept 

Gci[ 

1460.437 

-408.6899 

Gay 

1439.269 

-403.156 

Ni 

1908.567 

-504.769 

Ny 

1883.116 

-500.607 
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Interstitials per Ion per cm^ (Int/Ion/cm^) 



Figure 3.15. Linear Fitting for 2.0-MeV Protons in AlGaN 


Table 3.7. Interstitial and Vacancy Distribution Fitting Functions in the AlGaN 
Region of the Device 


Lines for 2.0-MeV energy levels. 


Atom 

Slope 

T-intercept 

Alj 

5015.562 

-2143.11 

Aly 

4723.538 

-2015.279 

Gai 

18310.142 

-7809.917 

Gay 

18157.719 

-7742.635 

Ni 

31059.048 

-13205.667 

Ny 

30677.914 

-13042.54 
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CHAPTER 4: 
Model Implementation 


4.1 Impurity Scattering Mobility Implementation 

The final piece of the mobility term developed in this work was for the elfect of impurity 
scattering. As stated, our device was not intentionally doped, so an impurity scattering 
term that is dependent on concentration of dopants returns a value that has no elfect on the 
initial model. Because of this, the tuning of the model to match the characteristics seen in 
the non-irradiated device was carried out without an impurity dependent term. 


To incorporate the elfects of impurity scattering in the absence of the first term of 
Albrecht’s mobility equation, the Brooks-Herring model. 


1 Qr/3r2 

l^BH ^imp 


(4.1) 


was used, where a and (3 are user defined scaling factors, Nimp is the impurity concentration 
in cm“^ and T is the temperature in K. The Brooks Herring model can be shown to provide a 
more accurate depiction of the impurity dependence of the scattering in a GaN 2DEG [25]; 
however, the true term should also depend upon the 2DEG density n 2 deg, which is not 
included for reasons of convergence. The a and /3 terms were tuned to be equal to 0.4 and 
1.0 so the model would match the measured mobility percent change values in [23]. The 
result of this fitting is shown in Eigure 4.1. 
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Figure 4.1. Percent Change in ijL 2 deg versus Fluence 


The fit for the mobility is almost exact, while the simulation shows an over estimate 
for all the other properties this study reviewed. 

4.2 Model Evaluation Versus 2.0-MeV Experimental Re¬ 
sults 

The data on the change in characteristics of the HEMT irradiated with 2.0-MeV protons 
were performed in [26]. In [26] sheet resistance Rsh, 2DEG concentration n 2 DEG, electron 
mobility /i 2 DEG> on state resistance Ron, transconductance Gni,sat, threshold voltage Vth were 
measured versus fluence. The experimental results are shown in Figure 4.2. 
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Radiation Dose (cm' ) 


Figure 4.2. Percent Change Values after 2.0-MeV Proton Irradiation. Source: 
[26]. 


We matched the fluence levels in [26] for the 2.0-MeV radiation simulation. We com¬ 
pared simulation results to four of the characteristics listed in the key of Figure 4.2: 2DEG 
electron mobility /i 2 DEG> threshold voltage Vth, on state resistance Ron, and transconduc¬ 
tance Gm,sat- The results of the comparison between //2 deg is shown in Figure 4.1. The 
comparison between the model outputs and experimental results are laid out in Figures 4.3, 
4.4, and 4.5, respectively. 
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Figure 4.3. Comparison Between Model Output and Experimental Results 
for Pereent Change in Vth versus Fluenee 



Figure 4.4. Comparison Between Model Output and Experimental Results 
for Pereent Change in Ron versus Fluenee 
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Figure 4.5. Comparison Between Model Output and Experimental Results 
for Pereent Change in Gni,sat versus Fluenee 


From Figures 4.3, 4.4, and 4.5, as stated above, the model overestimates the pereent 
change in comparison to the shifts seen in the physical device. This could be because of the 
annealing out defects in the physical device. The physical devices were not tested during 
the radiation event, while the model essentially freezes defects in place and does not allow 
annealing out. We expect the trap profile to begin to decrease over time as the vacancies and 
interstitials anneal, reducing the resultant damage. The model does not account for this. 

The IV characteristic output by the model after simulated 2.0-MeV proton irradiation 
to a fluenee of 6.0 x 10^"^ cm“2 are shown in Figures 4.6, 4.7 and 4.8 along with the 
un-irradiated model outputs. These plots clearly show the threshold voltage shift in the 
device towards zero at both 100.0 mV and 10.0-V drain voltage values. This matches the 
results seen in physically radiated devices, discussed further in Section 4.3. We also see the 
current saturate at a much lower value in the Id versus Vd curve shown in Figure 4.8, which 
is something that can easily affect the performance of a component utilizing one of these 
devices. 


37 








Drain Current (A) 



Figure 4.6. Vg versus 1^ Comparison, Vd = 100.0 mV 



Figure 4.7. Vg versus Id Comparison, Vd = 10.0 V 
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Figure 4.8. Id versus Comparison, Vg = -3.0, -2.0, -1.0 and 0.0 V 


4.3 Trap Occupancy and Modulation of the Electric Field 

In order to determine the physieal origin of the observed shifts in the deviee eharaeteristies, 
the internal strueture of the deviee at varying fluenee was examined. A review of the strue- 
ture files produeed by the simulation revealed ionized donor and aeeeptor trap oeeupaney 
mismateh in both the GaN and AlGaN regions of the deviee. This ean be seen in Figures 
4.9 and 4.10 for the AlGaN region and Figures 4.11 and 4.12 for the GaN region of the 
deviee. Both aeeeptor trap defeets maintained higher ionization levels than their donor 
eounterparts. 
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Figure 4.9. Trap and Donor Occupancy Levels (cm in the AlGaN Barrier 
Layer at Fluence level of 10^^cm“^ 



Figure 4.10. Trap and Donor Occupancy Levels (cm in the AlGaN Barrier 
Layer at Fluence level of 6.0 x lO^'^cm”^ 
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Figure 4.11. Trap and Donor Occupancy Levels (cm in the GaN Barrier 
Layer at Fluence level of 10^^cm“^ 



Microns 


Figure 4.12. Trap and Donor Occupancy Levels (cm in the GaN Barrier 
Layer at Fluence level of 6.0 x lO^'^cm”^ 


41 






















The trap and donor ionization difference creates a build up of negative charge in the 
AlGaN and GaN regions of the device. This charge buildup screens out the generation of 
the 2DEG, which results in the threshold voltage shift towards zero. After radiation, less 
charge is needed to turn the device off as the 2DEG is not as concentrated, resulting in 
decreased current for a given bias condition and decreased transconductance. 

The second effect is seen in the spread of the electric field, as shown in Eigures 4.13 
and 4.14. The electric field under the gate on the drain side of the device spreads out over a 
larger area at higher lluence levels. This increase in width decreases the field and increases 
the breakdown voltage of the device. The effect on the field under the drain edge of the gate 
can be seen more quantitatively in Eigure 4.15, where the field under the drain side of the 
gate is shown. 



109.6 110 110.4 

Microns 


Eigure 4.13. Electric Eield Distribution at a Eluence of 10^^ cm ^ 
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Figure 4.14. Electric Field Distribution at a Fluence of 6.0 x lO^'^cm ^ 



Microns 


Figure 4.15. Electric Field Distribution at Fluences of 10^^ cm ^ and 6.0 x 
lO^^crn-^ 
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4.4 Model Evaluation at Additional Energy Levels 

Using the additional radiation profiles derived from SRIM, we predieted the effeets that 
might be seen in an aetual radiation event at 5.0, 10.0, 20.0 and 40.0 MeV to similar fluence 
levels. The eomparison between the pereent ehange in key deviee eharaeteristies and proton 
energy levels ean be seen in Figures 4.16, 4.17, 4.18 and 4.19. The idea that the higher the 
energy level of the proton the more fluence is required to inflict a similar amount of damage 
to the device was expected to hold; however, the increase in fluence needed to shift device 
characteristics by a similar amount decreases as particle energy level increases. This is 
unexpected, as decreasing particle interaction should make it harder to damage the device. 
As the particle energy increases, we expect to see fewer interactions with lattice atoms, 
but when interactions occur, the results will be more extreme than at lower energy levels. 
Again, this is an issue of transferred energy. A 40.0-MeV proton has a lot more energy to 
transfer to a PKA than a 5.0-MeV proton. Where a 5.0-MeV proton may knock three atoms 
from their original location, a 40.0-MeV proton can knock out 50. Also, if an event occurs 
with a higher energy level proton, the proton now has less energy and chances of interacting 
with other particles increase. 



Figure 4.16. Percent Change in ^i 2 deg versus Fluence for Varying Proton 
Energies 
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Figure 4.17. Percent Change in Vth versus Fluence for Varying Proton Energies 



Figure 4.18. Percent Change in Ron versus Fluence for Varying Proton Ener¬ 
gies 


45 




















Figure 4.19. Percent Change in gm,sat versus Fluence for Varying Proton 
Energies 
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CHAPTER 5: 
Conclusions 


In this study, we have successfully created a TCAD model that utilizes Monte Carlo sim¬ 
ulated proton irradiation damage levels for several different proton energy levels while 
incorporating a customized calculation for the electron mobility in the 2DEG to simulate 
high-dose proton irradiation in a GaN-on-Si HEMT. 

Generally, the model over-predicted the damage that was observed in previous NFS 
and NRE experimental work. We did not over-predict the mobility shift in radiated devices 
because the mobility was used to tune the coefficients in the customized mobility calculation. 
One issue with using this metric to tune the model is that the devices that were physically 
tested were not immediately characterized after irradiation. This allows the device time 
to anneal out defects that formed during the irradiation. This can result in a shift back to 
less damage in the device. The model effectively freezes the defects in place in the lattice 
without any function to allow for them to anneal; however, the model provides a good 
baseline model for the exploration of radiation damage in GaN-on-Si HEMTs. 

Continued research into GaN-on-Si devices is vital to the Department of Defense (DoD) 
as more and more hardware systems adopt these HEMTs. They strengthen our equipment, 
but the most important part of strength is knowing weaknesses. This research increases our 
to understanding of these powerful devices that are now being utilized in active systems. 

5.0.1 Future Work 

The Silvaco software does not have specialized formulas which depend upon 2DEG mobility 
such as those proposed in research by Eisesivdin [25]. Eor future work on this project, more 
physically exact terms for mobility should be incorporated, such as the effects of roughness 
scattering. Eor the HEMTs that were used in this study, prior work by Greenlee [32] 
documented the visible shift in atoms at the boundary layer between the AlGaN and GaN 
layers, which has an effect on the 2DEG in that region through roughness scattering. 
Additional research should also examine various ways to model the electron mobility shift 
caused by PCE without physically altering the model. More research is also needed in 
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this simulation to model the diffusion of defects throughout the lattice in order to identify 
mechanisms for void sinks causing defects like the voiding seen under the gate of the 
device by Wade [10]. Lastly, work by Pearton identified different defect types created 
by different energy levels of proton radiation. Additional research should be devoted to 
identifying these trap levels and incorporating them into models for specific energy levels of 
proton irradiation. This model can serve as a valuable tool for predicting the effects in the 
experiments that should follow, 40.0-MeV proton irradiation of NRL GaN-on-Si ffEMTs 
to fluence levels greater than 10^^cm“^. This tool will also prove valuable for upcoming 
neutron irradiation testing. 
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APPENDIX 


A.l SRIM Simulation Values 

Table A.l. Interstitial and Vaeaney Fitting Funetions in the GaN Region of 
the Deviee 


Atom 

Slope {{defects/cm)/{ion! urn)) 

T-intereept 

Gcij 

395.816 

-83.213 

Gay 

648.499 

-198.918 

N, 

402.650 

-35.729 

Ny 

841.037 

-244.639 


5.0-MeV Protons 


Table A.2. Interstitial and Vaeaney Fitting Funetions in the AlGaN Region of 
Deviee 


Atom 

Slope {{de f ects/cm)/{ion//urn)) 

T-intereept 

All 

918.307 

-388.069 

Aly 

102.064 

-19.034 

Guj 

4299.410 

-1816.631 

Gay 

2817.0904 

-1117.765 

Ni 

10080.524 

-4296.504 

Ny 

2837.467 

-1042.749 


5.0-MeV Protons 


49 






Table A.3. Interstitial and Vacancy Fitting Functions in the GaN Region of 
the Device 


Atom 

Slope {{de f ects/cm)/{ion/ jam)) 

T-intercept 

Gcij 

315.347 

-91.402 

Gay 

310.406 

-89.853 

Ni 

437.557 

-123.811 

Ny 

432.481 

-123.114 


10.0-MeV Protons 


Table A.4. Interstitial and Vacancy Fitting Functions in the AlGaN Region of 
Device 


Atom 

Slope {{de f ects/cm)/{ion/[im)) 

T-intercept 

Al[ 

A13.033 

-193.481 

Aly 

453.97 

-185.344 

Gaj 

2358.752 

-983.149 

Gay 

2034.923 

-840.816 

Ni 

3345.905 

-1365.073 

Ny 

2752.402 

-1105.875 


10.0-MeV Protons 


Table A.5. Interstitial and Vacancy Fitting Functions in the GaN Region of 
the Device 


Atom 

Slope {{de f ects/cm)/{ion/[im)) 

T-intercept 

Gaj 

0 

31.988 

Gay 

0 

31.968 

Ni 

0 

45.575 

Ny 

0 

45.085 


20.0-MeV Protons 
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Table A.6. Interstitial and Vacancy Fitting Functions in the AlGaN Region of 
Device 


Atom 

Slope {{de f ects/cm)/{ion/ jam)) 

T-intercept 

All 

276.451 

-112.589 

Aly 

400.732 

-166.478 

Gcij 

1258.994 

-520.669 

Gay 

1550.988 

-647.579 

Ni 

949.152 

-352.061 

Ny 

1454.471 

-572.491 


20.0-MeV Protons 


Table A.7. Interstitial and Vacancy Fitting Functions in the GaN Region of 
the Device 


Atom 

Slope {{defects/cm)/{ion//urn)) 

T-intercept 

Gaj 

0 

17.561 

Gay 

0 

17.238 

Ni 

0 

25.968 

Ny 

0 

25.489 


40.0-MeV Protons 


Table A. 8. Interstitial and Vacancy Fitting Functions in the AlGaN Region of 
Device 


Atom 

Slope {{de f ects/cm)/{ion/nm)) 

T-intercept 

All 

-79.039 

40.018 

Aly 

-128.416 

61.496 

Gai 

681.631 

-282.472 

Gay 

686.527 

-284.880 

Ni 

790.314 

-310.940 

Ny 

844.579 

-335.516 


40.0-MeV Protons 
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A.2 SILVACO Simulation Values 


Table A.9. Albrect Model Values. Adapted from [33] 



0 

Bn 

5 

Cn 

0.4 

No 

lOi" 

To 

300.0 

Ti 

1065.0 


Values A, B and C were calculated based on //„ of 1265 from [34] with A constantly 
zero. 


A.3 Trap Levels in GaN 


Table A. 10. Defects Introduced by Various Types of Irradiation in GaN. 
Adapted from [20]. 


Material 

type, dopant 

concentration 

(cm“3) 

Radiation 

type 

Carrier 

removal 

rate, 

(cm“l) 

Defect levels 

(activation en¬ 
ergy (eV) from 

conduction 

C or Valence 

V), defect 

type (donor D, 
acceptor A) 

Defect pro¬ 
duction rate 

(cm“l) 

Proposed 

identity 

n-GaN, lO^^ 

Protons, 

260 

0.13, C, ? 

30 

ERl 


2MeV 


0.16, C, ? 

400 

ER2 




0.2, C, A 

600 

ER2 

n-GaN, lO^^ 

Protons, 

100 

0.2, C, ? 


ER3 


150 keV 


0.25, C, ? 


Complex of ER3 




0.32, C, ? 


Complex of ER3 




0.45, C, ? 


Complex of ER3 




0.6, C, A 






0.8, C,D 


Gai 
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Table A.l 1. Stated Energy Levels for Defects in GaN. Adapted from [20]. 


Atom 

Vacancy or Interstitial 

Acceptor or Donor 

Deep or Shallow 

Energy Level 

Ga 

Interstitial 

Donor 

Deep 

Ec - 0.8 eV 

Ga 

Vacancy 

Acceptor 

Deep 

E, + 1.0 eV 

N 

Interstitial 

Acceptor 

Deep 

Ec-l.O eV 

N 

Vacancy 

Donor 

Shallow 

Ec - 0.06 eV 


A.4 Code 

A.4.1 Silvaco Deckbuild Code 

go atlas simflags="-P 12" 

set radenergy = 2 

set fluence = 6el4 

set talgan = 17.5e-3 

set tgan = Q.8 

set talgan2 = 0.517+Q.258 

set tain = Q.427 

set tsi = 35® 

set tni = Q.Q3 

set taug = Q.3 

set tausd =0.33 

set tsin = Q.133 

set tair = 0.366 

set Isin = .133 

set Igate = 7 
set Is = 100 
set Isg = 3 
set Igd = 13 
set Id = 100 

set Tern = 21 
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set pz = 7.2el2 


mesh width=lQ® 

x.mesh location=0 spacing=$ls/5 

x.mesh location=$ls spacing=$ls/10® 

x.mesh location=$ls+®.5*$lsg spacing=$lsg/l® 

x.mesh location=$ls+$lsg-®.1 spacing=$lgate/18 

x.mesh location=$ls+$lsg spacing=$lgate/7® 

x.mesh location=$ls+$lsg+®.5*$lgate spacing=$lgate/l® 

x.mesh location=$ls+$lsg+$lgate spacing=$lgate/7QQ 

x.mesh location=$ls+$lsg+$lgate+®.1 spacing=®.®2 

x.mesh location=$ls+$lsg+$lgate+®.5*$lgd spacing=$lgd/l® 

x.mesh location=$ls+$lsg+$lgate+$lgd spacing=$ld/l®® 

x. mesh location=$ls+$lsg+$lgate+$lgd+$ld spacing=$ld/5 

y. mesh location=-l’'($tni+$taug+$tsin) spacing=$tsin/2 
y.mesh location=-l*($tausd) spacing=$tausd/l® 
y.mesh location=-l*($tni) spacing=$tni/l® 

y.mesh location=® spacing=$talgan/l® 

y.mesh location=®.99*$talgan spacing=$talgan/l®® 

y.mesh location=®.999*$talgan spacing=$talgan/l®®® 

y.mesh location=$talgan spacing=$talgan/l®® 

y.mesh location=$talgan+®.QQ2*$tgan spacing=$tgan/25® 

y.mesh location=$talgan+®.®2*$tgan spacing=$tgan/2®® 

y.mesh location=$talgan+®.2*$tgan spacing=$tgan/2Q 

y.mesh location=$talgan+$tgan spacing=$tgan/6 

y.mesh location=$talgan+$tgan+$talgan2 spacing=$talgan2/5 

y.mesh location=$talgan+$tgan+$talgan2+$taln spacing=l 

y.mesh location=$talgan+$tgan+$talgan2+$taln+$tsi spacing=$tsi/4 

region number^l material=AlGaN x.mole=®.27 x.min^® x.max=$ls+$lsg+$lgate+$lgd+$ld 
y.min=® y.max^Stalgan 

region number=2 material=GaN substrate x.min^® x.max=$ls+$lsg y.min=$talgan 
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y.max=$talgan+$tgan 

region number^S material=AlGaN x.mole=0.7 y.min=$talgan+$tgan 
y.max=$talgan+$tgan+$talgan2 

region number=4 material=AlN x.min^® x.max=$ls+$lsg+$lgate+$lgd+$ld 
y.inin=$talgan+$tgan+$talgan2 y.max=$talgan+$tgan+$talgan2+$taln 
region number^S material=Si x.min=® x.max=$ls+$lsg+$lgate+$lgd+$ld 
y.inin=$talgan+$tgan+$talgan2+$taln y.max=$talgan+$tgan+$talgan2+$taln+$tsi 

region number^G material=gold x.min^® x.max^Sls y.min=-l*$tausd y.max=® 
region number^? material=nitride x.min=$ls x.max=$ls+$lsg y.min=-l*$tsin y.max=® 
region number^? material=nitride x.min=$ls+$lsg-®.1 x.max=$ls+$lsg 
y.inin=-1*($tni+$taug+$tsin) y.max=-l*$tsin 

region number^? material=nitride x.min=$ls+$lsg+$lgate x.max=$ls+$lsg+$lgate+®.1 
y.min=-1*($tni+$taug+$tsin) y.max=-l*$tsin 

region number^? material=nitride x.min=$ls+$lsg x.max=$ls+$lsg+$lgate 
y.inin=-1*($tni+$taug+$tsin) y.max=-1*($tni+$taug) 
region number^? material=nitride x.min=$ls+$lsg+$lgate 
X.max=$ls+$lsg+$lgate+$lgd y.min=-l*$tsin y.max=® 

region number^S material=nickel x.min=$ls+$lsg x.max=$ls+$lsg+$lgate 
y.min=-l*$tni y.max^® 

region number=9 material=gold x.min=$ls+$lsg x.max=$ls+$lsg+$lgate 
y.min=-1*($tni+$taug) y.max=-1-'Stni 

region number^ll material=gold x.inin=$ls+$lsg+$lgate+$lgd 
X.max=$ls+$lsg+$lgate+$lgd+$ld y.min=-l*$tausd y.max=® 
region number=12 material=GaN x.inin=$ls+$lsg x.max=$ls+$lsg+$lgate 
y .inin=$talgan y .max=$talgan+$tgan 

region number^lS material=GaN x.min=$ls+$lsg+$lgate 
X.max=$ls+$lsg+$lgate+$lgd+$ld y.min^Stalgan y.max=$talgan+$tgan 
region number=14 material=air x.min^Sls x.max=$ls+$lsg-®.1 
y.min=-1*($tni+$taug+$tsin) y.max=-l*$tsin 

region number=14 material=air x.min^® x.max^Sls y.min=-l*($tsin+$tausd) 
y.max=-l*($tausd) 

region number^lS material=air x.inin=$ls+$lsg+$lgate+0.1 
X.max=$ls+$lsg+$lgate+$lgd y.min=-l*C$tni+$taug+$tsin) y.max=-l*$tsin 
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region number^lS material=air x.inin=$ls+$lsg+$lgate+$lgd 
X.max=$ls+$lsg+$lgate+$lgd+$ld y.min=-1*($tsin+$tausd) y.max=-1*($tausd) 


electrode name^gate x.min=$ls+$lsg x.max=$ls+$lsg+$lgate y.min=-l*$tausd y.max=® 
electrode name^source x.min=® x.max^Sls y.min=-l*$tausd y.max=® 
electrode name^drain x.min=$ls+$lsg+$lgate+$lgd x.max=$ls+$lsg+$lgate+$lgd+$ld 
y.inin=-l*$tausd y.max^® 
electrode name^substrate bottom 

doping n.type conc=lel8 uniform x.min=® x.max=$ls y.min=® y.max=$talgan 
doping n.type conc=lel8 uniform x.min=$ls+$lsg+$lgate+$lgd 
X.max=$ls+$lsg+$lgate+$lgd+$ld y.min^® y.max^Stalgan 

doping f.doping="GaN_N_I.lib" trap acceptor e.level^l.® degen.fac = 2 
sign=le-15 sigp=le-15 

doping f.doping="GaN_N_V.lib" trap donor e.level=3.34 degen.fac = 2 
sign=le-15 sigp=le-15 

doping f.doping="GaN_Ga_V.lib" trap acceptor e.level=2.4 degen.fac = 2 
sign=le-15 sigp=le-15 

doping f.doping="GaN_Ga_I.lib" trap donor e.level=2.6 degen.fac = 2 
sign=le-15 sigp=le-15 

doping f.doping="AlGaN_N_I.lib" trap acceptor e.level=l.® degen.fac = 2 
sign=le-15 sigp=le-15 

doping f.doping="AlGaN_N_V.lib" trap donor e.level=3.34 degen.fac = 2 
sign=le-15 sigp=le-15 

doping f.doping="AlGaN_Ga_V.lib" trap acceptor e.level=2.4 degen.fac = 2 
sign=le-15 sigp=le-15 

doping f.doping="AlGaN_Ga_I.lib" trap donor e.level=2.6 degen.fac = 2 
sign=le-15 sigp=le-15 

#doping f.doping="AlGaN_Al_V.lib" trap acceptor? e.level=2.4? degen.fac = 2 
sign=le-15 sigp=le-15 region^S 

#doping f.doping="AlGaN_Al_I.lib" trap donor? e.level=2.6? degen.fac = 2 
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sign=le-15 sigp=lG-15 region^S 


#Values taken from 8-3, 8-4 from pg 534 of atlas users manual using GaN/2 and 
AlN/2 values from CH6 of the users manual. Exponent term was not halved, 
material region=l affinity=3.6 tcon.power tc.npow = .317 tc.const = 1.381 
material region=2 affinity=3.932 tcon.power tc.npow = .28 tc.const = .65 
material region=3 affinity=2.78 tcon.power tc.const = 1.5535 tc.npow = .688 
material region=4 affinity=2.03 tcon.power tc.npow = 1.64 tc.const = 1.425 
material region=5 affinity=4.05 tcon.power tc.const = 1.48 
material region=6 tcon.power tc.const = 2.0 
material region=8 tcon.power tc.const = 1.58 

material region=12 affinity=3.932 tcon.power tc.npow = .28 tc.const = .65 
material region=13 affinity=3.932 tcon.power tc.npow = .28 tc.const = .65 

contact name^gate work=4.6 surf.rec 
contact name=source work=3.6 resistance=10 
contact name=drain work=3.6 resistance=10 

models print k.p fermi srh 
models lat.temp heat.full 

interface charge=$pz s.s y.min=$talgan y.max^Stalgan 

thermcontact number=l elec.number=4 temperature=273.15+$Tem 

mobility region=2 f.tofimun="MobO.lib" 
mobility region=12 f.tofimun="Mobl.lib" 
mobility region=13 f.tofimun="Mob2.lib" 

output con.band val.band charge int.charge polar.charge e.mobility 
method block newton itlimit=25 maxtraps=2Q nblockit^lOQ carriers=2 


57 



A.4.2 Tofimun Mobility Code 


GaN Region 0 

#include <stdio.h> 
#include <stdlib.h> 
#include <math.h> 
#include <ctype.h> 
#include <malloc.h> 
#include <string.h> 
#include <templatG.h> 


* ATLAS Parser Function Template 

* ATLAS Version 5.21. l.C 

* c 1993 - 2015 SILVACO Inc. 

* All rights reserved. 


V 

/-. 

* General field dependent mobility model for electrons. 

* Total Field Mobility (parallel and perpendicular field components) 

* Statement: MATERIAL/MOBILITY 

* Parameter: F.TOFIMUN 

* Arguments: 

* Eperp [in] - perpendicular electric field (V/cm) 

* Na [in] - acceptor concentration (/cm^3) 

* Nd [in] - donor concentration (/cm^3) 

* nconc [in] - electron concentration (/cm^3) 

* Eparl [in] - parallel electric field (V/cm) 

* TL [in] - lattice temperature (K) 

* xcomp [in] - x-species fraction (0-1) 

* ycomp [in] - y-species fraction (0-1) 
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* *mun [return] - hole mobility Ccm^2/Vs) 

* *dmundep [return] - derivative of *mun wrt Eperp 

* *dmundEparl [return] - derivative of ’-mun wrt Eparl 

* *dmundl [return] - derivative of *raun wrt TL 

* *dmundn [return] - derivative of *mun wrt nconc 


V 


int tofimunCdouble Eperp,double Na,double Nd,double nconc,double Eparl, 
double TL,double xcomp,double ycomp,double *mun,double *dmundep, 
double *dmundepar,double *dmundl,double *dmundn) 

{ 

double A=0.25*2.61e-4, B=4*2.9e-4, C=4*170e-4, N®n=lel7, T®n=30®.®, tln=l®65.®, 
alphas®.4,beta = 1.®,Ecrit=2®e3 ,N1=2.® ,Vsat=3e7; 
double invmu,DD®,DD1,DD2,invmuALB,h,g,hPrime,gPrime,MuField; 
double N_imp=Na+Nd; 

invmuALB=A’'(N_imp/N®n)*pow(TL/T®n,-1.5)+B*pow((TL/T®n),1.5)+C*pow((exp(tln/TL)-l),-1); 
DD®=-1.5''A*(N_imp/N®n)*pow(TL/T®n, - 2.5); 

DD1=B *(1.5)"pow((TL/T®n),®.5)"(1/T®n); 

DD2=C*tln''powC(exp(tln/TL)-l) , -2)’'(pow(TL, -2)) ; 

invmu=invmuALB; 

h=l+pow((Eparl/Ecrit),N1); 

hPrime=CNl/powCEcrit,Nl))’'powCEparl, (Nl-1)) ; 
g=(l/invmu)+Vsat''(powCEparl, CNl-l))/pow(Ecrit ,N1)) ; 
gPrime=(Vsat/pow(Ecrit,N1))*(Nl-l)*pow(Eparl,(Nl-2)); 

MuField=g/h; 

*mun = MuField; 

*dmundepar = ®; 

//(gPrime/h)-(1/invmu)*(hPrime/h); 

*dmundep=®.®; 


59 




-'dmundl = (DDO+DD1+DD2) / (invmu*invmu) /h; 
*dmundn = 0/(invmu*invmu)/h; 

returnC®); 


GaN Region 1 

Code is the same as seen in GaN Region 0 with the constant values changed. Only the 
constant values have been listed here. 

double A=0.25*2.61e-4, B=1.35*2.9e-4, C=1.35*170e-4, N0n=lel7, TOn^SOQ.Q, tln=l®65.0, 
alphas®.4,beta = 1.®, Ecrit=2®e3, Nl=2, Vsat=3e7; 


GaN Region 2 

Code is the same as seen in GaN Region 0 with the constant values changed. Only the 
constant values have been listed here. 

double A=®.25*2.61e-4, B=2*2.9e-4, C=2*17®e-4, N®n=lel7, T®n=3®®.®, tln=l®65.®, 
alphas®.4,beta = 1.®,Ecrit=2®e3 ,N1=2.® ,Vsat=3e7; 


A.4.3 Vacancy / Interstitial Formation Code 

AlGaN A1 Interstitial Code 

#include <stdio.h> 

#include <stdlib.h> 

#include <math.h> 

#include <ctype.h> 

#include <malloc.h> 

#include <string.h> 

#include <template.h> 
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* ATLAS Parser Function Template 

* ATLAS Version 5.21. l.C 

* c 1993 - 2Q15 SILVACO Inc. 

* All rights reserved. 


V 

/-. 

* Position dependent net doping. 

* Statement: DOPING 

* Parameter: F.DOPING 

* Arguments: 

* X location x (microns) 

* y location y (microns) 

* ’-nnet net doping concentration (per cc) 

V 

int doping(double x,double y,double *nnet) 

{ 

double talgan = 0.0175; 
double tgan = 0.8; 
if(y<0){ 

*nnet = 0; 

} 

else if(y>=0 && y<=talgan){ 
double flu=8e+14; 

*nnet= flu*5015.561905*(y+0.43)-2143.1099108; 

} 

else{ 

*nnet=0; 

} 

return(O); 

} 
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